Low density symmetry energy effects in the neutron star crust properties 
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Abstract 

The form of the nuclear symmetry energy E s around saturation point density leads to a different crust-core 
transition point in the neutron star and affect the crust properties. We show that the knowledge about E s 
close to the saturation point is not sufficient, because the very low density behaviour is relevant. We also 
claim that crust properties are strongly influenced by the very high density behaviour of E s , so in order 
to conclude about the form of low density part of the symmetry energy one must isolate properly the high 
density part. 
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1. Introduction 

One of the most intriguing quantity in the de- 
scription of nuclear matter in neutron stars is the 
symmetry energy E s , which is defined as follows 

E nuc (n,x) =V{n)+E s {n)a 2 + 0{a A ) (f) 

where E nuc (n, a) represents the energy of nucleonic 
matter per baryon as a function of baryon num- 
ber density n and the isospin asymmetry a, where 
a = (n n —n p )/n and the neutron and 

proton densities. At the saturation point density 
no = O.f 6 fm -3 the value of symmetry energy corre- 
sponds to the ct4 parameter in the Bethe-Weizsacker 
mass formula, and takes the value E s (jiq) = 30 ± 
I MeV. Isoscalar part of interactions is represented 
by the isoscalar potential V(n) which is mainly re- 
sponsible for the stiffness of the Equation of State 
(EoS). 

Density dependence of E s is however highly un- 
certain both below and above saturation point Uq. 
This dependence is one of the goals of the exper- 
imental investigations carried on the radioactive 



beam colliders 0, Q ■ This kind of facilities allow for 
research of nuclear matter with large isospin asym- 
metry. The analysis shown in [4, 3, 5] put some con- 
straints on the slope and curvature of E s (n) around 
no however we are still far from the final conclusion 
about the global shape of the symmetry energy. 

The role played by the E s in the context of vari- 
ous neutron star observables was emphasized in Q 
and more detailed analysis was made in Q . One of 
the first approach to the crust-core transition was 
presented in § where the stability considerations 
were performed. This kind of analysis with an im- 
proved nuclear model was later used in [9j. In this 
work authors suggested the different critical den- 
sity in different nuclear models comes from their 
discrepancy in the neutron matter description, e.g. 
in the symmetry energy form. The direct connec- 
tion of E s to the crust-core transition point was 
shown in [10( and possible phase separation in the 
inner core was analysed in In this work we 

would like to go along this line to emphasize that 
the very low density behaviour of the symmetry en- 
ergy is essential for the crust-core transition point. 
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It is especially interesting taking into account the 
recent experimental results 12|, |13[ which show the 
symmetry energy still takes large values at densi- 
ties very much below uq. This result is in contrast 
to common conviction coming from various theo- 
retical approaches that states the E s goes almost 
linearly to zero for low densities. So, the conse- 
quences of the symmetry energy with such unusual 
feature seems to be worth to be seen. 

2. The crust-core transition point 

In order to estimate the crust-core transition 
point the bulk instability conditions were applied. 
It is based on the vanishing compressibility which 
signals the one-phase system is unstable against the 
density fluctuations and must split into two phases. 
It was shown in [l(| that the compressibility un- 
der constant charge chemical potential - is the 
proper quantity in case of matter under beta equi- 
librium. It is defined as follows 



dP 

dn 



(2) 



(i 



where P is the total pressure (nucleons + leptons). 
When nuclear matter contribution to the total en- 
ergy is described by Eq. (fTJ) then the stability con- 
dition takes the form 



K u , = 



n 2 (E'Ja 2 + V 
2a 2 E' 2 n 2 



2n(E[a 2 + V') 



E s 
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(3) 



where the role played by the symmetry energy E s is 
apparent. For densities typical for the neutron star 
core the compressibility is positive but decreases 
as density decreases. For some density n c , located 
below no, the compressibility vanishes. It means 
the charge fluctuations are not stable - below n c the 
matter cannot exist in one phase and must separate 
into two phases which may form a crystal lattice. 
However, the answer, where is situated the inner 
edge of the crust, is much more complicated. 

There is no unique method to determine at what 
density the solid crust starts to form. Here we 
want to discuss this issue. The critical density for 
vanishing compressibility n c {K^) represents the ab- 
solute limit for the homogeneous and neutral sys- 
tem. It means the phase splitting must occur be- 
fore reaching the n c . Indeed, the energy for the 
two-phase system becomes smaller than for the one- 
phase system before the density approaches the 



critical value n c (X /i ) [3, Q. The boundary for 
two phase coexistence, n c (l O 2) is then located 
slightly above n c (K^) as it is shown in Fig|T] Both 



crust lattice 



pasta 
phases 



=0.1 n 



liquid core 



K >0 

M 



=1 n„ 



n (K) 



n(l~2) 



Figure 1: Schematic picture of the different estimation for 
the crust-core transition in the neutron star. The gray re- 
gion corresponds to < 0, the uncertain hatched range 
indicates region where pasta phases may occur. 

n c (K^) and n c (l ■<-> 2) are quantities derived in 
the bulk approximation, i.e. without the finite size 
effects included. The finite size structures require 
inclusion of the additional forms of energy like the 
Coulomb and surface energy. Competition between 
them, leads to structures with various dimension- 
ality (rods, plates) usually called pasta phases j22[. 
The presence of these exotic structures is however 
model dependent, moreover mechanical properties 
of rods and plates resemble rather liquid crystals 
than solid state [23|. The shear tensions in all 
spatial directions can be supported only by a 3- 
dimensional lattice of 0-dimensional structures, e.g. 
the lattice of almost point-like nuclei immersed in 
the less dense medium, e.i. the gas of electrons and 
dripped neutrons. So, for the most reliable position 
of the crust edge one should take the point where 
the nuclei start to deform into very long structures 
or, keeping their spherical form, dissolve in the sur- 
rounding medium. 

The finite size effects may also be included in the 
stability consideration of one-phase system as was 
done in The most important ingredient in this 
analysis appeared to be the Coulomb energy. Be- 
ing repulsive, the Coulomb force stabilizes fluctua- 
tions of charge and makes matter more resistive to 
formation of charged clusters. The corresponding 
critical density n c (Q) is then moved to lower den- 
sities, always below n c (if M ) (see FigfTJ. For some 
kind of interactions (e.g. Skyrme-like in [l8j]) the 
n c (Q) almost coincides with the point where nuclei 
disappear, so it may be treated as the very likely 
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Figure 2: The shape of symmetry energy for different models 
in the low and high density region (the inset). 



estimation of the crust-core transition. 

In this work we stay with the simplest condition, 
n c (K^) what means we get upper bound on the 
crust mass and thickness. The difference between 
n c (K^) and n c (Q) is not so large, usually smaller 
than 10%. The comparison between these quanti- 
ties was already done in Q and for modern nuclear 
models in [19ll2l| . The condition © corresponds to 
the cross-point of the border of spinodal region at 
zero temperature with beta equilibrium condition 
in stellar matter as it was considered in 1201. 



3. Nuclear models 

Different models of nuclear interactions, while be- 
ing consistent at saturation point uq = 0.16 fm -3 , 
leads to large discrepancies for EoS at higher densi- 
ties. In order to estimate the role played by E s we 
used a model which have the same isoscalar part of 
interactions V[n) and different form of the isovector 
part E s (n) The V[n) taken from Prakash et.al. fl5| 
correspond the compressibility at saturation point 
K = 240 MeV. 

Various shapes of symmetry energy, where taken 
from MDI model Q and named a, b, c, d as in [Icj |. 
The symmetry energy in MDI family presents large 
discrepancies both in the low and in the high den- 
sity region. Also the slope E' s (no) and curvature 
E" (no) at the saturation point take different val- 
ues. The only common property is the value of the 
symmetry energy at n y . Especially, the large differ- 
ences at high densities lead to a differences in the 
global stellar parameters like its radius and com- 
pactness. Those global parameters indirectly affect 
the crust thickness and its contribution to the total 



Table 1: The E B (n) parametrization for k models: 0,5,10. 
The last columns presents the critical density. 

a p 7 E s (0) [MeV] n c 

247.08 40.93 48.836 0.080 

181.98 103.43 15.502 5 0.073 

116.87 165.93 -17.831 10 0.065 



star mass what obscure the picture. In order to iso- 
late the low density properties of E s from the rest 
we use a specific parametrization which presents the 
same shape of E s at saturation point and very tiny 
differences at high density. It has the following form 



E s = 



au 2 + f3u + 7 V" + E s (0) 
u 2 + 1 



(4) 



where u — n/rioo. The 71^ control the behavior at 
very high density and was taken here = 1 fm -3 
to get intermediate values of E s between extremal 
cases (a and d) in the MDI family. The remain- 
ing parameters a, j3, 7 were derived to ensure the 
same slope L — 80 MeV and the stiffness K asy — 
—450 MeV at no while E s at n — > approaches 
three different values 0, 5 and 10 MeV. The shape of 
symmetry energy for all models are shown in Figj^l 
Above models of nuclear interaction cannot be 
use to obtain EoS for very low densities, in the crust 
region. In order to obtain the EoS int the full range, 
up to the edge of the star, they have to be completed 
with the EoS for the crust itself. Here we used well 
established results from [l6| and 17 1. 



4. Results 

In this section the results are presented. In the 
FigO the mass-radius relation is shown. One may 
observe huge differences in the MDI family. These 
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Figure 3: The total star mass versus star radius. 
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Figure 4: The star crust thickness versus the total star mass. 



differences are caused only by the symmetry energy 
shape because the isoscalar part is the same for all 
EoSs. The symmetry energies in the MDI family 
diverge enormously at higher densities. There are 
no results for the a model because its E s takes so 
large negative values at high density that makes the 
pressure bounded from above and leads to unstable 
EoS. The rest, b-d models, present large differences 
in the compactness of neutron star giving differ- 
ence in the surface gravitation. Hence any crustal 
parameters like thickness, mass etc. are changed 
not only by the form of E s at the crust region but 
also by its values in deep core region. For the k 
family E s is almost the same at high density and 
the M-R relations for fcO-fclO models are similar to 
each other. The differences between them slightly 
increase with the dropping mass what is natural as 
the crust occupies more an more volume of the star 
and differences in the EoS at low density are more 
pronounced. The effect of the high density E s form 
is apparent when one compare critical densities for 
different models. In the MDI family, for b, c, d one 
get n c = 0.092,0.095,0.160 (see [13]) whereas for 
fc0,fc5,fcl0: n c = 0.080,0.073,0.065 (see TabH|). 
The difference in critical density for b and c is less 
then in the k family but the difference in the crust 
thickness (see FigSJ) is much larger than one could 
expect from the n c position. Surely it comes from 
the different compactness of the star as a whole. In 
the c case the star is more compact and its crust 
is more squeezed by the gravity. The effect is even 
better seen for the crustal fraction of the moment 
of inertia, FigfS] As one can note for the typi- 
cal 1.5M star the relative differences among b, c, d 
models reach hundreds of percents whereas for the 
k family are not so large. At the end of the dis- 
cussion we focus on the k models themselves. The 



Figure 5: The crustal fraction moment of inertia versus the 
total star mass. The lower bound, 1.4%, for I cr ust/I coming 
from observations is indicated by the horizontal line. 



symmetry energy for this family was constructed 
to ensure the same behaviour at high density and 
keep the same slope and curvature at the saturation 
point in agreement with recent results from experi- 
mental results. The only difference was the asymp- 
totic value of the E s at vanishing density n 0. 
We have probed values 0, 5 and 10 MeV. Consid- 
ering the shapes of -E^in the Figf2] one may con- 
clude that the E s (0) is very essential quantity for 
the crust properties. The symmetry energies for 
the k family are almost overlapping in the no re- 
gion and below and their discrepancies seems to be 
negligible in comparison to b-d lines, however they 
effect in quite large differences for the crust thick- 
ness and moment of inertia. We also observe that 
higher £^(0) systematically makes the crust thin- 
ner and its contribution to total moment of inertia 
takes very low values. 

In order to refer our results to observations we in- 
dicated the upper bound on the crustal moment of 
inertia for the Vela pulsar coming from analysis per- 
formed by Link in 24J. The line for fclO model goes 
below 1.4% in very wide range of stellar masses. 
We do not know the Vela pulsar mass, so in that 
sense the model cannot be verified, but if we trust 
experimental results that E s (0) « 10 MeV it would 
mean the Vela pulsar mass should have its mass not 
greater than 1M . 



5. Summary 

We perform the analysis on how the very low 
symmetry energy behavior effects the basic neutron 
star crust parameters. Although the crust encom- 
passes the densities below no, also the high den- 
sity behavior of the symmetry energy influences its 



4 



properties. In order to be able to conclude about 
the low density part we constructed a special fam- 
ily of symmetry energies. They presented the same 
shape at saturation point and above but differ when 
density goes to zero. It was shown that unusual 
property like the non-vanishing symmetry energy 
at zero density leads to a very thin crust with very 
small contribution to the total moment of inertia. 
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